(^ 

0^ 

^ ; DIFFUSING RECONSTITUTING DIMERS: A SIMPLE 

^ ■ MODEL OF BROKEN ERGODICITY AND AGEING 

(N 

r— ni DEEPAK DHAR 

Q . Theoretical Physics Group, Tata Institute of Fundamental Research, 

^ I Homi Bhabha Road, Mumbai 400005 (India) 

_. - We consider a model of assisted diffusion of hard-core particles on a 

^_> I line. In this model, a particle can jump left or right by two steps to 

. . an unoccupied site, but only if the the site in between is occupied. 

jrt I We show that this model is strongly non-ergodic, and the phase space 

breaks up into mutually disconnected sectors. The number of such 
sectors increases as the exponential of the size of the system. Within 
' O . a sector, the model can be shown to be equivalent to the Heisenberg 

5 ! model, and is exactly soluble. We study how time-dependent autocor- 

relation functions in the model depend on the sector. We also show 
how this strong breaking of ergodicity ( which can be thought of as an 
infinity of conservation laws) affects the hydrodynamical description 
of the long wavelength density fluctuations in the system. Finally, we 
^vj ' study the effect of allowing transitions between sectors at a very slow 

Q^ ' rate. 
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f-^ ■ 1 Slow and fast processes 






If a system is weakly coupled to a heat bath at a given "tempera- 
ture", if the coupling is indefinite and not known precisely, if the 
coupling has been on for a long time, and if all the "fast" things 
'^ I have happened and all the "slow" things not, the system is said to 

be in thermal equilibrium. R. P. Feynmanu 
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The Boltzmann-Gibbs method of calculating properties of a sys- 
^ . tern of large number of interacting particles by defining a suitable 

ensemble of systems, and replacing the calculation of time-averages 
of obsevables by their phase space averages, is very powerful tech- 
nique in statistical mechanics. To calculate the equilibrium proper- 
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ties of a system of a large number of degrees of freedom, all we have 
to do is to calculate the corresponding partition function, and thus 
the free energy. All quantities of interest can then be determined 
by taking derivatives of this free energy with respect to suitable 
conjugate fields. Of course, this calculation is quite nontrivial, 
involves evaluation of many-dimensional integrals, and can be ex- 
actly done only for a very small class of soluble models. However, 
standard approximation techniques, such as weak- or strong- cou- 
pling perturbation expansions, or variational schemes can be used 
to get answers which are as close to the correct answer as desired 
away from critical points, and near critical points other schemes 
such as renormalization group, e— expansion etc. can be used to 
determine the behavior of different quantities. Thus, it is perhaps 
not an exaggeration to say that the general principles of determin- 
ing equilibrium properties of systems are quite well understood by 
now. 

However, many systems of interest, such as diamonds or ordi- 
nary glass, are not in thermal equilibrium, even though their prop- 
erties do not seem to change significantly with time for the time- 
scale of experiments. These systems are said to be in metastable 
states. It is well realized that these states have much higher free 
energy than some other stable state of the system, and if the sys- 
tem could easily explore all parts of the phase space available to 
it, such states would occur with negligible probability. They owe 
their existence to the fact that the rate of transition to the stable 
phase can very low, and the system is not ergodic. 

A precise microscopic specification of all positions and momenta 
of all the atoms, to describe the state of a small piece of glass is 
impossible. Clearly, we can, at best, hope to give a probabilistic 
description of the state. However, the ensemble of microscopic 
states having finite probability must to restricted to those which are 
'nearby', and are accessible to the system within times comparable 
to the experimental time-scales. Such an ensemble in these systems 



is much smaller than the corresponding microcanonical ensemble. 
It is easy to see that phase space volume of accesible states of 
such a system occupies a negligible fraction of the full phase space, 
and this fraction decreases as exp{—aV), as the volume of system 
V increases. The coefficient a seems to have only a very weak 
dependence on the time scale used to define the set of 'accessible' 
states. Of course, a must become zero as this time tends to the 
Poincare recurrence time of the system. 

However, it has proved difficult to define this concept of re- 
stricted ensembles in a way which allows calculation of such parti- 
tion functions, and thus determine the average properties of, even 
model systems. This is what we shall try to do here for a very 
simple toy model called the diffusing reconstituting dimers (DRD) 
model. The ideas developed in this paper have been discussed in a 
recent paper □. In this present paper, our aim is to argue that this 
is a very nice and tractable, hence useful model for glassy dynam- 
ics, while the earlier paper dealt more with its integrability, and 
connection to other soluble models. We shall not attempt to review 
here the very large body of work which already exists dealing with 
the phenomena of ageing and relaxation in spin glassesQ. What is 
presented is a somewhat personal perspective. 

The key observation here is note that the processes operating 
in a glass may be divided into two classes: fast and slow. Fast 
processes are those whose time scales are much shorter than that 
of observation, e.g. the vibratory motion of atoms about their 
mean positions. Slow processes are those involving rearrangement 
of atomic configurations, particle- or vacancy diffusion, creep etc. 
To define a quasi-equilibrium state of glassy systems, we have to 
imagine that "all the fast things have happened", and the slow 
things have not even started. To be precise, we assume that the 
system dynamics is such that the slow processes do not occur at all. 
Then the dyanmics is nonergodic. The phase space breaks up into 
a large number of disconnected pieces, called sectors. The number 




Figure 1: A schematic representation of the decomposition of phase space into 
disconnected sectors. The points denote different configurations. Full lines 
show allowed fast transitions between sectors. Different sectors are enclosed in 
rectangular boxes. Dotted lines show slow transitions between different sectors. 



of such sectors increases exponentially with the size of the system. 
We shall call such system many-sector decomposable (MSD), and 
the ensembles which sum over only one of these sectors as pico- 
canonical ensembles. 

For the toy model we consider, partition functions correspond- 
ing to such pico-canonical ensembles can be explicitly computed. 
These differ from sector to sector. As the number of sectors is of 
the order of exp{V), a precise specification of the sector is too mi- 
croscopic a detail, and not relevant. In an experiment, one can, 
at best, hope to give some general probabilistic characterization 
of the sector. This presumably would depend on the history, and 
preparation of the system, e.g. cooling schedules etc.. In a the- 
oretical calculation, this implies that once the free energy in the 
pico-canonical ensemble is determined, it must then be further av- 
eraged over sectors with a suitable weight for each sector. 

Finally, we shall show how we can describe the phenomena such 



as ageing in glasses by taking into account the slow processes, as 
providing a stochastic dynamics in which the system jumps from 
one sector to another in time, and thus brings back ergodicity ( or 
tries to ) at very large times. 



2 The Model 

The model we shall choose to address these issues is a very simple 
one. We consider N hard-core point particles on a linear chain 
of L sites. There are no interactions amongst the particles except 
the hard-core interaction. We assume that the system evolves in 
continuous time time by a local Markovian dynamics. We shall 
choose the transition rates to satisfy the detailed balance condition, 
so that the rate from a transition, say from configuration C to C, 
is the same as the rate of the reverse transition C" to C. Then, 
in the steady state all accessible configurations occur with equal 
probability. 

Now for a more detailed specification of the allowed transitions. 
We assume that the particles can diffuse to nearby sites, but only if 
there is another particle nearby: a case of assisted diffusion. More 
precisely, a particle at site i can jump left or right by two steps 
to the sites i±2, if that site is unoccupied, and the in-between 
site i±l is occupied. If allowed, this transition occurs at rate 1. 
We represent the configuration of the system by a binary string 
01100100.... of length L, where the i*'^ bit gives the occupation 
number Ui of the site i. Then this process can be represented by a 
'chemical' equation 

110^ >011 

An alternate way of viewing the system is to say that isolated 
particles cannot diffuse. But two particles at adjacent sites can 
diffuse together one step to the left or right with rate 1. These 
pairs of particles will be called dimers. However, these pairings are 



not for ever, and the dimers can 'reconstitute'. Thus, for example, 
in the sequence of transitions 

..11010.. ^..01110.. ^..01011.. 

the middle particle is first paired with particle on the left, and then 
in the second transition with the particle on the right. We are thus 
considering a system of Diffusing Reconstituting Dimers (DRD). 

A third possible way of looking at the model is to think of the 
diffusion as that of O's by two steps left or right if the intervening 
site is not a zero. This is then not assisted diffusion, but too many 
zeroes will 'hinder' each other's diffusion. 

The assumption in the model monomers have a much lower dif- 
fusivity than dimers is not easily realized in real physical systems. 
However, instances where this happens are not unknown: for ex- 
ample, platinum dimers on some surfaces have higher mobility than 
monomers. 

3 Many Sector Decomposition and Conservation Law of 
the Irreducible String 

For of linear chain of L sites, the phase space consists of 2^ distinct 
configurations. It is however clear that the dynamical rules are such 
that the system is not ergodic. For example, the as the diffusion 
process conserves the number of particles, a system starting with N 
particles at some time will never be found in a configuration with 
a different number of particles. The phase space thus breaks up 
into mutually disconnected subsets, called sectors. The stochastic 
evolution of the system takes it from one configuration to another, 
within the same sector, but no transitions are allowed between 
different sectors. 

The conservation of particle number implies that the total phase 
space can be broken into (L + 1) disconnected subsets, each sub- 
set specified by a different value of the total particle number. Such 
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conservation laws are well known in equilibrium statistical mechan- 
ics, and there effect is easily taken into account by defining suitably 
restricted ensemble : say the canonical ensemble in which the num- 
ber of particles is fixed. In addition, usually, one can prove that 
there is an equivalence of different ensembles ( except at special 
points like when the system is at a phase transition point). So one 
can even work with ensemble having a variable number of parti- 
cles, so long as the chemical potential is chosen so that the average 
density comes out right. 

But this is not all. It is easy to see that if we break the linear 
chain into two sublattices (even and odd), then number of particles 
on each sublattice is separately conserved. Thus the phase space 
has to be broken up further. Most of the (L -|- 1) subspaces have to 
further subdivided into roughly L parts, so that the total number of 
disconnected regions of phase space is now of order L^ . This would 
imply that to describe the long-time steady state of the system 
correctly, one would need not one, but two chemical potentials, 
one for each sublattice, and adjust them so that the initial state 
density on each sublattice is correctly reproduced. 

So long as the number of conserved quantities is finite, one can 
continue this process, and eventually end up with a system de- 
scribed by conventional Gibbs measure on a suitably defined 'mi- 
crocanonical' ensemble in the phase-space, or equivalently a 'canon- 
ical' ensemble where the conservation laws are ignored, but their 
effects taken into account by including a finite number of thermo- 
dynamic potentials^. If there are m conserved quantities, and each 
takes order L distinct values, the phase space has order U^ distinct 
sectors. 

However, it is easy to see that in the DRD model, the breaking 
of ergodicity is much more severe. In fact the number of discon- 
nected sectors increases as exponential of L. This is easy to see. 
Consider a configuration in which there are no adjacent sites both 



occupied such as 

...0010100010010000101. 



From the rules of the DRD dynamics, such a configuration can- 
not evolve at all, and remains the same at all times. The number 
of such configurations is easily shown to grow exponentially with 
L. [It is easy to show that this grows as yU^, where /i is the golden 
mean =1.618...] 

^ 111010010 

^ 101110010 
> 101011010 

^ 101001110 

101001011 

Figure 2: A non-trivial sector with L=9, and five states within the sector. 
The arrows indicate allowed transitions between these. For this sector the 
irreducible string is 1010010. 

Consider now a less trivial case of single diffusing dimer in a 
frozen background. This is shown in fig. 2. We see that while the 
movable I's change in time, as the dimer reconstitutes as it moves, 
there is exactly one movable dimer in each such configuration if you 
want to move to the right, and exactly one (possibly different) if 
you want a leftward move. After an order L moves in one direction, 
one reaches the end, and no further moves are possible. So there 
are about L configurations in each such sector. The number of 
sectors again grows as jj,^. 



It is possible to continue this process, and construct sectors in 
which 2, 3, .. movable dimers are present. To get a full description 
of the phase-space decomposition, we need to develop a systematic 
method for doing this. This is most conveniently done in terms of 
a construction called the irreducible string (IS), which we describe 
next. 

With each of the 2^ possible configurations of the DRD model, 
we attach a binary string, called the IS corresponding to that con- 
figuration, constructed as follows: We take the L-bit binary string 
of Hi specifying the configuration. We read this string from left to 
right until the first pair of adjacent I's is encountered. This pair 
is deleted, reducing the length of the string by 2. We repeat this 
process until no further deletions are possible. The resulting string 
is the IS for the configuration. For example, for the binary string 
01001110110, the irreducible string is 0100100. 

By construction, for each configuration, there is a unique IS. 
However, many different configurations may yield the same IS. 
The usefulness of this construction stems from the following ob- 
servation: Two different configurations belong to the same sector, 
if and only if they both have the same IS. This is easy to see. 

Firstly, we note that we need not assume that the pair of I's 
to be deleted has to be the first occurrence of consecutive I's when 
read from left to right. We get the same IS, whatever the choice 
and order of pairs 11 deleted, so long as the final string has no pairs 
of I's in it. [For example, from the configuration 101110, we get 
the same IS 1010, whether the second and third I's, or the third 
and fourth I's are deleted.] Then, if we get configuration C from 
C in one elementary step, deleting the diffusing dimer first, we see 
that both C and C have the same IS. Thus, the IS is a constant 
of motion. 

Secondly, if two configurations C and C have the same IS, then 
there exists a sequence of allowed transitions of the DRD model, 
which takes C to C'.To prove this, notice that we can push any dif- 
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fusing dimer in C as far right as possible, and get a standard config- 
uration in the sector, in which the configuration is the IS followed 
by all the dimers. Thus, for example, the standard configuration 
corresponding to the configuration 101110 is 101011. Correspond- 
ing to the configuration 01001110110, it is 01001001111. If C and 
C have the same IS, they can both be transformed into the same 
standard configuration, and as the DRD rules are reversible, into 
each other. So they belong to the same sector. 

The conservation of the IS thus gives us a convenient way to 
label the different sectors of the DRD model. To each allowed 
IS corresponds a unique sector. In addition, this decomposition 
of phase space into disjoint sectors using the IS as a constant of 
motion is maximal in the sense that there are no other independent 
constants of motion. 

As a side remark, let us note that the number of independent 
constants of motion is not a well-defined concept for a discrete sys- 
tem. For a classical mechanical system with n-dimensional phase 
space, we say that there are r independent constants of motion if 
the trajectories are confined to (n-r) dimensional manifolds. But 
for systems with discrete phase space, the trajectory consists of a 
finite number of discrete points, and its dimension is undefined. 

In fact, for such a system, given two independent constants of 
motion /i and I2, one can always construct a third single constant 
of motion J3, such that both /i, and I2 are known functions of I^. 
This result was proved in the context of quantum mechanics, by von 
Neumann eI (for a recent discussion see WiegertH). His argument is 
quite simple, but does not seem to be as well known as it should 
be. Consider Ji, I2 as two given commuting quantum mechanical 
observables over a finite dimensional Hilbert space. Since /i and I2 
are mutually commuting, they can be simultaneously diagonalized. 
Then the quantum-mechanical Hilbert space of all possible states 
of the system can be represented as a direct sum of subspaces, each 
of which corresponds to a pair (Ai, A2), where Ai and A2 are eigen- 
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values of the operators Ji and I2. We order these subspaces in an 
arbitrary order, and define the operator Is to be block diagonal in 
this basis, with eigenvalue r for all eigenvectors in the r*'* subspace. 
Then, by construction. This operator commutes with both Ji and 
/2, and its eigenvlaue r tells us the sector, and hence both Ai and 
A2. Thus, the operators Ji and I2 are functions of I^. 

Using the von Neumann construction, we can take any set of 
commuting constants of motion Ji, J2, /s, ..., and collapse them into 
a single constant. In general, if one can construct an observable 
which takes the same value for all configurations within a sector, 
but distict values in distinct sectors, then such an observable is a 
constant of motion and there can not exist any other constants of 
motion independent of it. Conversely, a single constant of motion 
having 2^^ distinct eigenvalues can be thought of as n independent 
constants of motion, each having only two distinct eigenvalues. 

For the DRD model, this implies that the IS can be thought 
of as a single constant of motion, equivalent to a large number 
of independent constants of motion: If we count the number of I's 
between the i*^ and (i + 1)*'^ zero in the binary string corresponding 
to the configuration, with time this number can only change in 
multiples of 2, hence its parity is conserved. If one is uncomfortable 
with constants of motion whose values are character-strings, the IS 
can be thought of as a binary integer. 

There is another way to represent the IS as a constant of motion, 
which makes contact with the theory of integrable models. We 
consider two noncommuting matrices A(O) and A{1) ( these can in 
general be complex n x n martices, but may be taken to be 2 x 2 
real matrices for our purpose here) which satisfy the relation 

A{l)A{l)AiO) = A{0)A{1)A{1) (1) 

For any configuration C of the DRD model with occupation num- 

11 



bers Hi, we define a matrix Il{C) by 

h{c) = n A( 



(2) 



where the matrix product is ordered so that larger j values are 
more to the right. Then, from the eq(l) it follows that this matrix 
product does not change by an elementary evolution step of the 
DRD model. This implies that II will not change in time even as 
C changes, and thus is a contant of motion. Eq(l) does not com- 
pletely determine A(0) and ^4(1), as these forms only 4 equations 
for 8 real parameters making the marices A{Q) and A{1). Thus we 
can choose a one parameter family of matrices A(0, A) and A(l, A) 
depending on a real parameter A. A simple specific choice is 



A(0,A) 



and 



" ] 


L-A A 




1 -1 




" 1 " 






-1 





A(1,A) 
Now II becomes a 2X2 matrix of the form 

Il{C,\)- 



/21(A) 



/12(A) 
/22(A) 



(3) 



(4) 



(5) 



whose elements In, /12, /21 and I22 are polynomials in A of degree 
at most L. Since this matrix is conserved for all values of A, it 
follows that all the coefficients of these polynomials are constants 
of motion. These coefficients in terms can be written out explicitly 
as functions of rii. 

4 Equivalence to a Multi-Species Exclusion Process 

The algorithm to determine the IS also specifies the positions of the 
characters in the L-bit binary string that remain undeleted. As the 
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configuration evolves, these positions will change, but the relative 
order of the characters remains unchanged. The DRD model can 
thus be viewed as a process of diffusion of hard-core particles on a 
line. Let the IS be of length L Then the number of particles is L 
Let Xi{t) be the position of the i^^ character of the IS (counting 
from left). Each particle carries with it a 'spin label' which is 
or 1, and is unaffected by the dynamics. Thus, our model is a 
special case of the Hard- Core Random Walkers with Conserved 
Spin (HCRWCS) models which may be defined for all 1-d models 
where a conserved IS can be defined QoEM. 

If the i^^ character of the IS is a 1, it must be followed by a in 
the IS, and also on the chain, this character must be immediately 
after it. so that we have 

X,+i(t)=X,,(t) + l 

for all times t, and the i*^ and {i + 1)*'^ walkers always move to- 
gether. In this case it is better to think of these as a single walker, 
which occupies two consecutive sites on the chain. Then there are 
two kinds of walkers: walkers with spin label occupying a single 
site, and walkers with label 10 which occupy two sites. In addition 
there are the sites not occupied by the random walkers, which are 
occupied by dimers 11, which again occupy two sites. Let us call 
these particles of type S,C and A respectively. We now show that 
the DRD dynamics can be viewed as a stochastic dynamics of a 
system of 3 species A, B and C of particles on a line. In the lat- 
ter process ( to be called the Exclusion process (ExP)), each site 
of linear chain is occupied by a single particle, which may be of 
type A, B, or C. We set up a one-to-one correspondence between 
the configurations of the DRD model, and of the ExP as follows: 
Read the binary string of the DRD from left to right. If the next 
character of the DRD model is 0, another character in the DRD 
string is read. We write a single character Aor B for the ExP string, 
depending on if the last chacter read was 1 or 0. Conversely, the 
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binary string corresponding to a configuration ABBCAABCB... 
of tlie ExP is obtained by a simple substitution A — i> 11, i? — > 10, 
C^O. 

We can now write tlie stocliastic evolution rule in terms of the 
ExP configurations. These are quite simple: A particle of type A 
can exchange position with a particle of type S or C at an adjacent 
site with rate 1. Type B and C cannot exchange positions with each 
other. This model thus is a special case of the k-species exclusion 
process defined by Boldrighini et allHl. In the latter, there are k 
distinct species of particles (k=2,3,4..). and particles can exchange 
places with adjacent particles with a rate which depends on both 
the species of the particles beinf interchanged. In our case k=3. 
The k-species model has been studied recently Derrida et alB, and 
by Ferrari et alli^. 

The conservation of IS in this language corresponds to the sim- 
ple statement that as B and C type particles cannot exchange 
places, their relative order is unchanged in time. Consider a string 
composed of three characters A, B and C which specifies a config- 
uration of the ExP process. Then deleting all the occurrences of 
A's in the string, we are left only with a string composed of 5's 
and C's which is conserved by the dynamics. 



5 Calculation of the Number and Sizes of Sectors 

In terms of the ExP model, it is quite straightforward to write down 
the formulas for the number of distinct sectors for the DRD model, 
and also the number of configurations in each sector. Consider 
free boundary conditions for convenience. In a sector having N^ 
particles of type A, Nb particles of type B, and Nc particles of 
type C, the total number of sites is L = 2Na + 2Nb + Nq. The 
number of I's in this case is 2Na + Nb, and the number of O's is 
Nb + Nc. The length of the IS is 2Nb + Nc- 

The number of distinct sectors is the number of distinct ways of 
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writng the IS consisting of Nb -B's and Nc C's. Since all choices are 
allowed, this is {Nb+ Ncy. / (NbINcI) ■ The number of configuration 
within a sector with a specified IS, say BCCCBBCBC...., is the 
number of ways we can place A^^ ^'s between (A''^ + Nc) -B's and 
C's. So this number is 

Q{Na, Nb, Nc) = (N^ + Nb + iVc,)/[iV^!(iV5 + Nc)\] (6) 

This determines the picocanonical entropy in each sector. To cal- 
culate the total number of sectors, we have to sum the number of 
sectors over all allowed values of Nb and Nc- 

Multiplying Q{Na, Nb, Nc) by the number of sectors with A^^ 
B-type and Nc C-type particles, we get the the total number of 
configurations having Na, Nb and Nc particles of type A,B,C re- 
spectively. Denote this by 'E:{Na,Nb,Nc)- Then we get 

E{Na,Nb,Nc) = {Na + Nb + Nc)\/[Na\Nb\Nc\] (7) 

Extremizing this with respect to Na, Nb and Nc, subject to the 
constraint that 2Na+'2Nb+Nc = L, we see that the this is extrem- 
ized for Na = Nb = Nc/2 = L/6. Hence if we pick a configuration 
at random, its length is most likely to be approximately 2L/3, and 
it will have approximately L/6 diffusing dimers. The typical num- 
ber of configurations in this sector is k'" where k = 2^/^.3"^'^. The 
number such sectors is approximately k'^, where k' = 3-'^/^.2~^/^. [ 
We check that kk,' = 2, so that the total number of configurations 
increases as 2"^, as it must.] 

Alternatively, we can proceed directly as follows: Let -^^(0) and 
-F^(l) be the number of distinct strings of length i which end in a 
and a 1 respectively. Since in the IS, a 1 must be followed by 0, by 
a may be followed by either or 1, these F's satisfy the recursion 
relations 

F,+i(0)=F,(0) + F,(l) (8) 

F,+i(l)=F,(0). (9) 
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Given the starting values -Fi(O) = -Fi(l) = 1 these recursions are 
easily solved, and we see the -^^(l) is the £*^ Fibonacci number 
which grows as /i^, where /x = (v^ + l)/2. For a line of length L, 
the largest possible value of £ is L. A sector with r diffusing dimers 
corresponds to i = L — 2r. Hence the total number of sectors is 

[L/2] 

Nl= J2[FL-2r{0)+FL_2r{l)] (10) 

r=0 

where [L/2] is the largest integer not greater than L/2. This shows 
that the total number of sectors also increase as fi^ for large L. 

For the DRD model, the probability weights in the steady state 
correspond to the rather trivial Hamiltonian H=0, and so the cal- 
culation above does not seem to have much interesting structure, or 
variation from sector to sector. But it does provide us with a sim- 
ple system where one can completely characterize different sectors, 
and calculate phase-space averages within a specified sector. 

6 The Rate matrix as a Quantum Spin Hamiltonian 

The technique of writing down a quantum-spin hamiltonian corre- 
sponding to a Markov process for a discrete system is by now quite 
well-known ll^. For completeness, we recapitulate it here briefly. 
Let P{C, t) be the probability that a classical system undergoing 
Markovian evolution is in the configuration C at time t. These 
probabilities evolve in time according to the master equation 

^E^ = E \-^ic ^ c) Pic t) + w{c' -. c) Pic t)] (11) 

where the summation over C is over all possible configurations of 
the system and W{C ^ C) is the transition rate from configuration 
C to configuration C . 

We define a Hilbert space, spanned by basis vectors |C), which 
are in 1 to 1 correspondence with the configurations C of the system. 
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A configuration with probability weight P{C, t) is represented in 
this space by a vector 

\P{t))=Y.PiC,t)\C) (12) 

c 

The master equation can then be written as 

^ ^ -HIPW) (13) 

This equation can be viewed as an imaginary-time Schodinger equa- 
tion for the evolution of the state vector \P{t)) under the action of 
the quantum Hamiltonian H. 

The 2^ configurations of the DRD model on a line of length L 
can be put in one-to-one correspondence with the configurations 
of a spin-1/2 chain of L sites. At the site i, the spin variable is 
erf = 2raj — 1, and takes values -|-1 and —1 if rii is 1 or respectively. 
It is then straightforward to write H as the Hamiltonian of the spin 
chain in terms of the Pauli spin matrices cxj. We find that H has 
local 2-spin and 3-spin couplings, and is given by 

^ = -^E[^-i-^m-l](l + ^-) (14) 

i 

This Hamiltonian is quite similar in structure to an integrable spin 
model with three-spin couplings proposed and solved by Barie^Ej. 
The Bariev model has the Hamiltonian given by 



The sign of U can be changed by the transformation af -^ erf, 
af —>■ erf, and af — > —erf. We set \U\ = 1. Then the Bariev model 
differs from our model through the term 

H-Hs = H' = -^Y: (<i^r+i - l) (1 + ^n (16) 

i 
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It is easy to see that these terms commute with the IS operator 
/^(A). Thus /l(A) also commutes with Hb) and provides an infinity 
of constants of motion of the Bariev model (in the special case 
\U\ = l). 

It is useful to write this Hamiltonian in terms of fermion oper- 
ators cl and Ci defined by the Jordan- Wigner transformation 



exp 



exp 



i-l 



cr; 



cr, 



(17) 



In terms of these fermion operators, H can be written as 



where rij = cfci is the number operator at site i. The first two 
terms of this Hamiltonian represent assisted hopping over an occu- 
pied site, and the last term describes two- and three- body poten- 
tial interaction between nearby sites. Similar hopping terms are 
encountered in a model studied by Hirsch in the context of hole 
superconductivity^. 

7 Conservation laws of the Quantum Hamiltonian 

While the evolution of probabilities in a classical Markov process 
could be cast as the evolution of the wavefunction of a suitably 
defined quantum mechanical problem, the concept of constants of 
motion is quite different in these two cases, as I explain below. 

A classical observable O is a function which assigns a value, 
say in terms of a real number 0{C) to each possible configuration 
C of the system. We say that (9 is a constant of motion when the 
value of the observable 0{Ct), associated with the configuration of 
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the system Ct at any arbitrary time t, is independent of t. In a 
quantum mechanical formulation, observables are represented by 
matrices O, and the quantity of interest is the expectation value of 
the observable O at time t, given by the formula 

{d), = {i\d\p{t)), (19) 

where (1| is the row vector with all entries 1. Clearly, in this 
language, a classical observable O corresponds to a quantum me- 
chanical operator O, which is diagonal in the natural basis of the 
system, i.e. when the basis vectors are \C). 

We shall call a quantum mechanical observable O a constant 
of motion, if the matrix O commutes with the Hamiltonian H. In 
this case it is easy to verify that its expectation value evaluated 
with respect to the probability vector \P{t)) does not change in 
time. (Note that we are using a different prescription for evaluating 
expectation values than in standard quantum mechanics.) Clearly 
the class of all possible quantum mechanical constants of motion 
(all matrices O) is much larger than that of the classical constants 
of motion (diagonal in the configuration basis). 

We have already seen that the IS is a classical constant of mo- 
tion. Using its definition, (Eq. 4.5), we define the diagonal matrix 
hiX) by 

h{X) = l[A{n,,X), (20) 



where fij = crtcyj is now an operator, and not a c-number. It is 
easy to verify that H commutes with /l(A) , for all values of A. 
Thus /l(A) provides us with a one parameter family of constants 
of motion. In fact, as discussed earlier, there are no additional 
independent classical constants of motion in our model. 

As H commutes with /l(A), it has a block-diagonal structure, 
with each block corresponding to a sector of the phase space, and 
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to a particular IS. The task of diagonalizing H then reduces to that 
of diagonahzing it in each of the sectors separately. 

Consider any one of these sectors. Let the IS X for this sector, 
in the 3-species exclusion process notation, be a string of Nb -B's 
and Nc C's in some order, say X = BCBBCBCC . . . The number 
of diffusing dimers in this sector then is A^^ = {L — 2Nb — Nc)/2. A 
typical configuration in this sector is specified by a string of length 
(A^^ + A''^ + Nc) with Na A^s interspersed between the characters 
of the IS, e.g. BCABAABCBCAC . . .. 

We define a chain of spin- 1/2 quantum spins Tj, with i ranging 
from 1 to L^ = (A^^ + Nb + Nc). For each configuration in the sec- 
tor X, we define a corresponding configuration of the r-spin chain 
by the rule that r^ = +1 if the j^^ character in the string specifying 
the configuration is A. For r? = —1, the corresponding character 
can be either B or C. But this degeneracy is completely removed 
by using the known order of these elements in the irreducible string 
X. Thus there is a one-to-one correspondence between the config- 
urations of the DRD model in the sector X, and the configurations 
of the r-chain with exactly Na spins up. 

The action of the Hamiltonian H on the subspace of configu- 
rations in the sector X looks much simpler in terms of the r-spin 
variables. It is easily checked that the evolution in terms of the 
r-variables is that of a simple exclusion process. The quantum 
Hamiltonian for this process is the well known Heisenberg Hamil- 
tonian 

Lr-l 

HHe^s = - E (^^ ■ ^^+l " ^) ' (21) 

1=1 

Note that Hueis looks different for sectors with differing L^'s. We 
have defined the transformation between the a— and r— spin vari- 
ables by an algorithmic prescription. It seems very difficult to do 
so by an explicit formula. The transformation is very complicated, 
and non-local. 

We would now like to construct a one-parameter family of oper- 
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ators that commute with Hhcis- This is a well-known construction 
for the Heisenberg model til. We use periodic boundary conditions 
for convenience. Define 2x2 matrices Lj{^) whose elements are 
operators acting on the spin Tj, and /i is a parameter 

L,(/.)=fl + ^t^ /^"^ I (22) 

and define 

t(/x) = rrriL,(/i). (23) 

Then it can be shown that E3 for all /i, [i' 

■r(/.),r(/.ol=o. (24) 



Or^ Oa 



0, for all r, s and 



Writing In r(/i) = E^i /^^^ we get 

J2 = Hhcis- Thus the set of operators {J} constitute a set of 
quantum mechanical constants of motion for the Hamiltonian Hueis 
in each sector separately, and hence for H. 

The Hamiltonian H has still another additional infinite set of 
constants of motion. These are related to the existence of an addi- 
tional symmetry in the model. Clearly, replacing a B-type particle 
by a C-type particle does not affect the dynamics of the exclusion 
process. Thus, the full spectrum of eigenvalues of H in any two 
sectors of the DRD model with different irreducible strings, but 
having the same values of Na and Nb + Nc-, is exactly the same. 
Such a symmetry may be viewed as a local gauge symmetry be- 
tween the B and C 'colors'. Changing the character in the IS from 
B to C (or vice versa) changes the length of L by 1. This symmetry 
therefore relates two sectors of the DRD model with different sizes 
of the system. 

The simplest inter-sector operators which preserve the total 
length of the chain L are operators which interchange two char- 
acters of the IS. Let Kj be the operator which interchanges the j*^ 
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and (j + 1)* characters of the IS. Then clearly, we have 



H,K, 



Oior j = lto Lr-1. (25) 



This 'color' symmetry implies that the spectrum of H is the same 
in all the sectors with the same value of L and Nj^. 

We have thus shown that the DRD hamiltonian shows the ex- 
istence of three different classes of infinity of constants of motion: 
the classical-mechanical constants of motion which may be encoded 
in the conservation of IS, the off-diagonal constants of motion of 
the Heisenberg model which can be used to diagonalize the martix 
H within a sector, and the constants of motion K which are also 
off-diagonal, but relate eigenvalues and eigenvectors in different 
sectors. 

8 Time-dependent Correlation Functions 

The static correlation functions in the DRD model are quite trivial 
because of the simple form of the hamiltonian chosen. The ef- 
fect of ergodicity- breaking and constrained evolution is best seen 
in the time- dependent density autocorrrelation functions. These 
show interesting variations from one sector to another. Such varia- 
tions occur despite the fact that, in each sector, there is a mapping 
between the DRD model and the simple exclusion process whose 
dynamics is known to be governed by diffusion. As we will see be- 
low, this mapping leads to a correspondence between tagged hole 
correlation functions in the two problems, but the form of auto- 
correlation function decays depends on the IS, and can be quite 
different. 

Consider a particular sector with IS I. Let the number of dif- 
fusing A, B and C particles in the equivalent 3-species exclusion 
process be Na, Nb and Nc respectively. In the IS, let h{k) be the 
number of B's to the left of the k-th hole.( A hole is a particle of 
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type B or C.) Evidently, the function b{k) specifies the IS com- 
pletely. 

Different configurations in the sector are obtained from different 
distributions of A's in the background of B's and C's. If there are 
Ofc A's to the left of the k'th hole, the location of this hole in the 
ExP and DRD problems is 

yExp{k,t) = k + ak{t), (26) 

yDRDik,t) = k + h{k) + 2ak{t). (27) 

Notice that h{k) (unlike a^) is time- independent. Defining a tagged- 
hole correlation function for the ExP and DRD problems in analogy 
with the conventional tagged-particle correlation functioiEa through 

a\t)=<[y{k,t)-y{kM^>^ (28) 

we see that the simple relationship 

ODRoit) = 2aExp{t) (29) 

holds. However, no simple, exact equivalence between the two 
problems can be established for tagged- particle correlations or single- 
site autocorrelation functions. This is because the transformation 
between spatial coordinates in the DRD and related ExP problems 
is nonlinear; a fixed site in the former problem corresponds to a 
site whose position changes with time in the latter. 

Let X and ^ denote spatial locations in the DRD and ExP prob- 
lems respectively and let a^ be the number of A's to the left of ^. 
Evidently, the number of holes (B's and C's) to the left of this site 
is (^ — a^). In the equivalent configuration in the DRD problem, 
each A corresponds to two particles, while the (^ — a^) B's and 
C's occupy a length 6(^ — a^) + E, — a^ which depends on the IS in 
question. Thus 

x = a^ + h{i- a^) + e, (30) 
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while the number of particles a^ to the left of x is given by 

Oa; = 2a^ + 6(^ — a^). (31) 

The transformation between the integrated particle densities a^ and 
a^ is therefore quite complicated. It depends on the IS through the 
function h and is highly non-linear. 

The correlation functions involving a^ are quite simple. The 
density-density correlation function for the ExP in steady state is 
defined as 

CE.p{i, t) =< n'(eo + e, io + t)n\U to) > ■ (32) 

where n' is a particle occupation number and an average over ^o 
and to is implicit. CexP satisfies the simple diffusion equation. 

^^^f^=V|C..p(e,t), (33) 

where V? is the discrete second- difference operator. For large t, 
therefore, CexP decays as t~^/^. Since a^ is a space integral over 
n'(^), correlation functions involving a^'s can be obtained as well. 
However, the change of variables from a^ to a^ (Eqs. 8.5 and 8.6) 
is difficult to perform explicitly. Nevertheless, we can determine 
the asymptotic behaviour of correlation functions CDRD{x,t), de- 
fined analogously to Eq. 8.7. Our technique for determining the 
long-time behavior of these correlations has been used earlier in 
discussing other models of the HCRWCS type &tlil3E2l. It relies on 
the fact that the slowest modes in the problem are diffusive motion 
of characters of IS, and the persistence of local density fluctua- 
tions of I's determines the decay of autocorrelation function at 
large times. We illustrate the technique below for various sectors. 
All these conclusions have been verifled by extensive Monte Carlo 
simulations!. 

The simplest sector is the one characterized by the IS 101010 .. . 
of length £ which is a finite fraction of the total length L. In this 
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sector, all the odd sites are always occupied and the dynamics on 
the even sites is that of the simple exclusion process. It is known 
that in the latter, the density- density autocorrelation function has 
a r~^/^ tail for large times r. Thus Cdrd{xo,t) is zero for xq on 
the odd sublattice, while it decays as r^^/^ on the even sublattice. 

Next consider the case where the IS consists of all zeros 0000 . . ., 
whose length is a nonzero fraction of L. Again, the t~^/^ decay of 
CexP then implies a similar behaviour for Cord- 

Now consider the general case of an IS whose j^^ character is 
aj = 0, 1. Let us assume that at time to, the site Xq is occupied 
by a particular character of the IS. Between the time to and to + 1, 
let the net number of dimers which cross the point xq towards the 
right be m (a leftward crossing being counted as a contribution —1 
to m) . For large times t, it is known that the distribution of m is 
approximated by a Gaussian whose width increases as t^^^, i.e. 



Prob(m|t) ~ — — exp 

v27rAf 



2 



(34) 



where A^ increases as t^^^ for large £^. If m dimers move to the 
right, a site occupied by aj at time to will now be occupied by 
%+2m at time to + t. Hence the autocorrelation function C{t) is 
approximated by 

oo 

C{r) = Yl "j"i+2m-Prob(m|t), (35) 

m=—oo 

where TxJaJ^^ is the average value of the correlations of characters 
in the IS averaged over j. This can have different values for even 
and odd j's as a particular element of the IS always stays on one 
sublattice. Define 7 odd (m) = a7a~7^ averaged over odd/even 

even ' J J 

sites. We have 



even ^-^' even 



Codd (r) = V'7odd (m) Prob(m|t). (36) 

fitrein ' ^ even \ 1 / \ / 
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If 7(m) is a rapidly decreasing function of m, then only small values 
of m contribute to C(r). Thus Cdrd varies as t~^^^ whenever 
correlations in the IS are short ranged. This will be so for most 
sectors of the DRD model. However, we can imagine preparing the 
system in way such that 'y{k), the Fourier transform of 7(m), varies 
as /c" , as |A;| — > with a a given adjustable parameter. In such a 
situation, the decay of Cdrd{t) with r will be power law with the 
exponent depending on a. 



9 Hydrodynamic Description of a Many-Sector Decom- 
posable System 

Usually, to describe the behavior of a system for large length or 
time scales, the hydrodynamical description is most appropriate. 
The number of fields that have to be used in such adescription 
depends on the number of conservation laws in the system. For 
example, usual hydrodynamics involves 5 densities ( of particle 
number, momentum and energy) corresponding to the fact these 
are the only locally conserved quantities. In systems with an in- 
finity of conservation laws, one would expect an infinity of such 
densities, which would make a hydrodynamical description impos- 
sible. Of course, well-known hydrodynamical equations such as 
the Korteweg de Vries equation are known to have an infinity of 
constants of motion. We would like to understand what is the re- 
lationship between existence of infinity of independent constants 
of motion and the hydrodynamical description. Our study of the 
DRD model does not allow us to answer this question fully, but 
suggests that the appearance of arbitrary functions in the hydro- 
dynamical equations of motion is a signal of partial integrability. 

In the DRD model, a hydrodynamic description would have to 
be in terms of coarse-grained density fields Podd{x, t) and peven{x, t) 
where a; is a continuous variable (0 < a; < L). These are the 
local densities of O's on odd and even sublattice sites respectively. 
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Equivalently, we can use the integrated density fields niodd and 
nieven respectively, where 

'modd{x,t) = podd{x',t)dx', (37) 

Jo 

and rrieven {x, t) is defined similarly. The equations of motion of rriodd 
and rrieven can be transformed into simple linear diffusion equations 
characterizing the exclusion process by a non-linear change of vari- 
ables from coordinates x to coordinates ^ (Eqs.30,31). 

The constraint that O's do not cross (equivalent to the restric- 
tion that the B and C particles of the 3-type exclusion process must 
maintain their relative ordering), implies that rrioddix, t) is a known 
function of meven{x,t) (and vice versa). Thus we can write 

modd{x,t) = f{meyen{x,t)), (38) 

where f{m) is a positive, non-decreasing function of its argument 
m, with /(O) = 0. Evidently / is determined completely by the 
irreducible string, as all O's reside on the IS. Alternatively, we may 
determine / completely if the initial occupations of each sublattice, 
Podd{x,t = 0) and Peven{x,t = 0) are known. The function /(m) 
which is specified by the initial conditions appears explicitly in the 
equations of motion as a constraint condition. This equation of 
constraint is a consequence of the conservation of the IS. 

In a more general context, we can imagine n coarse-grained 
density fields {(/)}, whose evolution equations are such that after 
suitable transformations of fields, we can define a field 0i, such 
that its evolution does not depend on the remaining fields, and 
hence can be written in the form 

^ = F(M. (39) 

while the other fields {0j}, J = 2 . . . n may involve couplings to 0i 
in addition to couplings to each other. If the evolution equation 
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for 01 is integrable, we can explicitly write (f)i{x,t) in terms of 
(f)i{x,t = 0). This would imply that the system has an infinity of 
conservation laws. In that case, the equations of motion of the other 
fields {4>j}, will involve an arbitrary function (f)i{x,t = 0). In the 
DRD model, the equation(38) is a constraint equation involving the 
hydrodynamical fields and an initial condition dependent function. 

10 Slowly Restoring Ergodicity 

We shall now try to describe the phenomena of ageing within the 
context of the DRD model. Thus, we would like to describe the 
system at time scales when the slow processes can not be fully 
ignored. These slow processes have rates which decrease sharply 
with temperature, and are negligible at temperatures lower than 
some 'glass temperature'. At larger temperatures, these processes 
would cause transitions between different sectors at finite rates, 
and restore ergodicity. 

For the DRD model, the simplest process that causes a tran- 
sition between different sectors ( and thus changes the IS), would 
be exchanging neighboring B- and C- type particles in the IS. This 
process is however not local in the DRD model, as in order to decide 
if a particular 1 at a site is part of an A- type or B- type particle of 
the ExP, we have to look to the left till a is encountered, and this 
may be many spaces away. The simplest process which changes the 
IS, and is local in the DRD model is 

...0010...^ ^...0100... (40) 

In terms of the ExP, this corresponds to the reversible transition 

...CCB...< >...CBC... (41) 

Let us assume that this process occurs at rate F, which is much 
less than 1. Then for times t << 1/F, we have broken ergodicity 
and the IS is a constant of motion. For times of order 1/F, the 
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system moves from one sector to another. However, in this case, 
the ergodicity is not completely restored. Of course, in this case 
Nb and A'',^ are constants of motion. But we are not concerned 
about O(L^) decomposition of phase space. If Nc > Nb, then it 
is easy to see that the (A''^ + Ncy./{NB\Ncl) sectors with having 
same value of A''^ and Nc get connected by the slow dynamics, 
and ergodicity is restored. If however Nb > Nc, then the C-type 
particles can only move about in a limited region of space, and not 
be able to go everywhere. 

For example, if in the IS , the number of -B's between two 
consecutive C's is everywhere at least two, it is easy to see that 
such IS's cannot evolve even under the slow dynamics. The number 
of such IS's is easily seen to grow exponentially with L. Thus this 
slow process merges sectors, but still an exponentially large number 
of them survive. The ergodicity is restored only partially. 

Consider an IS with a finite number of C's in a sea of 5's. Say 
the configuration is . . . BBBCCCBBBB .... Then it is easy to see 
that the leftmost C cant move at all, and the rightmost C by at 
most two characters to the right. In general, a configuration with 
n consecutive C's can spread out so that there are no adjacent C's, 
leaving the leftmost C fixed, and no further. 

It is straight forward to count the number of sectors in the 
presence of the slow process. For each of these big sectors, we can 
get a standard configuration by by pushing all C's as far to the right 
as possible. Then its structure is a random string formed with the 
substrings CB and B in arbitrary order. It is easily seen that this 
number is Nb\/[Nc\{Nb — A'c)!], which grows exponentially with 
the size of system. To study relaxation of auto-correlations in the 
presence of this slow relaxation is quite straight forward. For times 
t much less than l/F, the autocorrelation typically decays as t^^/*. 
However for larger times in case ergodicity is restored fully, it will 
crossover to a diffusive t~^/^ decay. In case Nb > Nc, the decay 
would continue to be t~^^^ for large times, but presumably with a 
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modified coefficient. Further study of tliis model is in progress. 



11 Generalizations and Relation to other models 

In summary, the DRD model is a very simple model for studying 
dynamics of MSD systems i.e. those systems showing strong break- 
ing of ergodicity, in which the phase space can be decomposed into 
an exponentially large number of mutually disconnected sectors. 
For this model, we can determine the sizes and numbers of these 
exactly. The sectors can be labelled by different values of a con- 
served quantity, the Irreducible String. The exact equivalence of 
the model to a model of diffusing hard core random walkers with 
conserved spin allowed us to determine the sector-dependent be- 
haviour of time-dependent correlation function in different sectors. 
In any given sector, we showed that the stochastic rate matrix 
was the same as that of the quantum Hamiltonian of a spin— 1/2 
Heisenberg chain (whose length depends on the sector), and thus 
demonstratated that it was exactly diagonalizable. 

The construction of the irreducible string in this model is very 
similar to other one- dimensional stochastic models with the many- 
sector-decomposability propertvstudied recently, i.e. the fc-mer 
deposition-evaporation modell^ciaEJ and the g-color dimer deposition- 
evaporation model BII2I. In all of these models, the long-time be- 
haviour of the time-dependent correlation functions can be deter- 
mined by noting that it is expected to behave qualitatively in the 
same way as that of the spin-spin autocorrelation function in the 
HCRWCS model with the same spin sequence as the IS in the cor- 
responding sector. 

However, the DRD model differs from earlier studied models 
in significant ways. In the trimer deposition evaporation model 
(TDE), the correspondence between the configurations on the line 
and the position of hard core walkers is many to one, unlike the 
present model where it is one-to-one. As a consequence, in the 
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steady state of the TDE model, all configurations of the random 
walkers are not equally likely, and one has to introduce an effective 
interaction potential between the walkers which is found to be of 
the form 

\/({X,}) = 5:/(X,+i-X,) (42) 

i 

where {Xi} axe the positions of the walkers, and f{x) increases as 
|lnX for large X. 

In the TDE model, the transition probabilities for the random 
walkers are also not completely independent of the spin-sequence of 
the walkers. In the g-color dimer deposition evaporation (qDDE) 
model, the color symmetry of the model implies that the dynamics 
of random walkers is completely independent of the spin sequence 
of the walkers, but the potential of interaction V is still present, 
which makes the problem difficult to study exactly. The DRD 
model is thus simpler than both the TDE and the qDDE models, 
and has the additional virtue of being exactly solvable in the sense 
that the stochastic matrix can be diagonalized completely. 

There are some straightforward but interesting generalizations 
of the model. Consider a general exclusion process with k types of 
particles. In this general model, if we assume that some types of 
particles cannot exchange positions (setting their exchange rate to 
zero), their relative order will be conserved and this can be coded in 
terms of the conservation of an IS. As a simple example, consider a 
model with 4 species of particles labelled A,B, C and D respectively 
with the allowed exchanges with equal rates 

AC ^ CA, AD ^ DA 

BC^CB BD^DB. ^^ 

This model again has the MSD property. It is easy to see that 
there are now two irreducible strings which are conserved by the 
dynamics. This is because the dynamics conserves the relative or- 
der of A and B type particles, and also of C and D type particles. 
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Figure 3: An example of a partially stuck configuration of the DRD model in 
two dimensions. The particles denoted by black circles are completely immo- 
bile, and cannot move under the DRD rules. The mobile particles ( denoted 
by white circles)can do so only along some horizontal or vertical lines. 
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In a string specifying the configuration formed of letters A,B,C and 
D, deleting all occurrences of A and B gives rise to an IS specify- 
ing the relative order between C and D type particles, which is a 
constant of motion. Similarly, deleting all occurrences of C and D 
characters, we get another independent IS which is also a constant 
of motion. In a specific sector, where both irreducible strings are 
known, the dynamics treats A and B particles as indistinguishable, 
as also C and D. Thus the dynamics is the same as that of the 
simple exclusion process (with only two species of particles), and 
is equivalent to the exactly solved Heisenberg model. 

The DRD model is interesting in higher dimensions also . For 
example, it is easy to see that on a square lattice in two dimen- 
sions, the number of totally jammed configurations increases expo- 
nentially with the number of sites in the system. All configurations 
with no two adjacent I's are totally jammed. These are just the 
configurations of the hard-square lattice gas modec3, whose num- 
ber is known to increase exponentially with the area of the system. 
One can also construct configurations in which almost all sites are 
jammed, except for a finite number of diffusing dimers, which move 
along a finite set of horizontal or vertical lines (see fig. 3). Thus the 
number of sectors in which only some of the particles are moving 
also increases as exponential of the volume of the system. How- 
ever, there is no equivalence to the 2-d Heisenberg model, and no 
obvious analog of the IS to label the disjoint sectors, and solving 
for the dynamics seems quite difficult. 

An interesting generalization of the model would be to the case 
when the hamiltonian is non-trivial, say an Ising model with nearest 
or nearest and next-nearest neighbor couplings. Then the allowed 
processes are assumed to be same as in the simple DRD model 
discussed above, but the rates are assumed to satisfy detailed bal- 
ance. In this case, in the steady state, realtive weights of differ- 
ent configurations are given by the Boltzmann weights, and one 
can determine the pico-canonical partition functions as functions 
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of parameters like temperature in any given sector. The relative 
weights of different sectors in the sector-averaging define another 
temperature parameter in the system. This may be identified as 
the temperature at which freezing occurs( quenching temperature). 
Details will be discussed in a future publication. 

Another generalization of the model would be to make the dif- 
fusion asymmetric. The corresponding 3-species exclusion process 
then becomes asymmetric, and belongs to the KPZ universality 
clascB. The invariance of dynamics under exchanging B and C 
type particles DD again reduces the problem to a simple asymmet- 
ric exclusion process of 2 sDecies, known to be exactly soluble by 
Bethe ansatz techniques c3'E3'El, and has a non-classical dynamical 
exponent 3/2. The correlation function of the asymmetric DRD 
model would map to somewhat complicated multispin correlation 
functions of the simple asymmetric exclusion process. How these 
would vary from sector to sector has not been studied so far. 

One intriguing feature of this model ( and other models with 
conserved IS) is the fact that we are able to get an infinity of clas- 
sical constants of motion by constructing a one parameter family 
of matrices which commute with each other and with the Hamilto- 
nian. A similar construction is used in integrable models, involving 
the R-matrix which is assumed to satisfy the Yang -Baxter equa- 
tion. The similarity of construction suggests that one can devise a 
more general R-matrix from which both the classical and quantum 
mechanical conserved quantities are deducible. Such a structure 
remains to be found. 
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